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The theory of hght scattering for a system of hnear molecules with anisotropic polarizabilities 
is considered. As a starting point for our theory, we express the result of a scattering experiment 
in VV and VH symmetry as dynamic correlation functions of tensorial densities pim{q) with I = Q 
and 1 = 2. I, m denote indices of spherical harmonics. To account for all observed hydrodynamic 
singularities, a generalization of the theory of Schilling and Scheidsteger ||^ for these correlation 
functions is presented, which is capable to describe the light scattering experiments from the liquid 
regime to the glassy state. As a microscopic theory it fulfills all sum rules contrary to previous 
phenomenological theories. We emphasize the importance of the helicity index m for the microscopic 
theory by showing, that only the existence of m = 1 components lead to the well known Rytov dip in 
liquids and to the appearance of transversal sound waves in VH symmetry in the deeply supercooled 
liquid and the glass. Exact expressions for the phenomenological frequency dependent rotation 
translation coupling coefficients of previous theories are derived. 
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I. INTRODUCTION 



In 1957 Rytov |^ predicted, based on a macroscopic phenomenological theory, that the depolarized light scattering 
should show a dip at zero frequency. Later in the 60 's it was indeed found experimentally |^ that the light scattering of 
a liquid in VH geometry, where the plane of polarization of the incident light and the scattered light are perpendicular 
to each other, shows a symmetric doublet at ± 1 GHz, with a minimum at zero frequency. However Rytov's theory 
only showed to be in superficial agreement with the experiments. 

In the following time a number of theoretical attempts have been made to understand the spectra. In general the 
intensity of the scattered light I^^ is proportional to the correlation function of the fluctuations in the components 
Q^S/ ((fj t) of the total polarizability tensor and to the incident light . Scattering in VH geometry is of particular 
interest since the direct (VV or HH) contributions are dominated by Brillouin peaks which correspond to propagating 
sound waves. In an ideal depolarized spectrum the Brillouin lines will be absent and therefore further information of 
the low frequency dynamics of the system can be obtained. 

It is most commonly believed that in molecular liquids, where anisotropic single molecule polarizability fluctuations 
are the main cause of depolarized spectra, the dip is caused by the coupling of molecular orientations to some part of 
the stress tensor. These theories can all be understood in the framework of the Mori-Zwanzig projection technique 
by using different variables for the projection scheme. Some primary variable is picked which couples via a memory 
kernel to other secondary variables. Until the end of the 70th however the generalized viscosities which occur in the 
memory kernel where assumed to be time independent. 

In 1969 Volterra 0| used the orientations as a primary variable to which the stress tensor couples which he believed 
to have - due to symmetry reasons - a nonzero matrix element with the total polarizability. Anderson and Pecora 
proposed in 1971 Q a theory which used only the symmetric part of the stress tensor as a secondary variable. Keyes 
and Kivelson in 1972 Q used the total polarizability as a primary variable and the momentum density as a secondary. 
Ailawadi et al. j7| in 1972 coupled the spin angular momentum density to the asymmetric part of the stress tensor. 

At that time it became evident that when reaching towards the supercooled regime additional features occurred 
e.g. additional site peaks or a central peak which could not be explained consistently by any of the above two variable 
theories. Therefore a number of theories were developed with additional variables (and an additional number of 
adjustable parameters). Quentrec in 1976 ||] used the whole second rank tensor of the orientations and Chappell 
et al. in 1981 Q coupled the momentum density and an unspecified symmetric tensor to the total polarizability. 
And also a four variable theory with 7 adjustable parameters was proposed ||l^. Wang in 1980 |l^ was the first to 
introduce time dependent generalized viscosities which can account for retardation effects in the memory functions. 



An approach which got refined and improved later on |12 13]. The authors of ||l3] and |l3| emphasize the importance 



of retardation - (memory -) effects in the coupling coefficients between the rotational and translational degrees of 
freedom. Since all the mentioned theories are phenomenological, it is difficult to decide, which one will give the 
"correct" description of the spectra. As an important result of our theory, which we are going to present, we will 
show, how the most complex of the mentioned phenomenological theories [O] can be derived microscopically as a 
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special approximation to a more complete set of equations. Thereby microscopic expressions for the phenomenological 
coupling constants in jl^] are presented. 

All of the above theories have in common, that they deal with the theory of the experimental technique (light 
scattering) and the dynamic of the system, which is measured in light scattering, on equal footing. One important 
deficiency of the early theories, is a too simplistic treatment for the dynamic of the liquid (e.g. assuming white noise 
spectra for the in general frequency dependent generalized viscosities |^), although the phenomenology for the light 
scattering part is correct. From a conceptual point of view, a clear distinction between these two topics should be 
made. The question which quantity is measured, can and should be answered first. In a second, and we would like to 
stress, independent step, a theory for the measured quantity can be formulated. This strategy was already pursued in 
the generalized hydrodynamic approach of fl^ , where - without specifying the scattering mechanism - a completely 
general treatment of the effect of hydrodynamic singularities and the influence of slow structural relaxations on light 
scattering spectra was achieved, by deriving formally exact expressions for frequency dependent Pockels constants 
and related Green Kubo relations. 

Also in our approach the two mentioned questions are clearly separated. But contrary to [p^, we will formulate a 
theory for a specific scattering mechanism, to get some more insight in the microscopic mechanisms. We will show in 
section^, that the formal structure of the theory is compatible with First we will derive the quantities, which 
are measured in an ideal light scattering experiments at linear molecules, where all interaction induced effects can be 
neglected. Then, we present a set of microscopic equations, which in principal enable us to calculate these quantities 
for a supercooled liquid close to the glass transition. In a final step we approximate our equations very drastically, 
keeping only the necessary ingredients for a qualitatively correct description of the light scattering experiments. This 
step is only for pedagogical reasons, to demonstrate the ability of the theory to reproduce light scattering spectra 
close to the glass transition. 

This paper is organized as follows: In section |l| we derive the direct contribution to the light scattering spectrum 
of a liquid of linear molecules. It will provide the starting point for the following theoretical considerations. 

In section [II we present a straight forward generalization of the molecular mode coupling theory ^ for linear 



molecules, which in addition to the orientational and translational degrees of freedom also describes the coupling 
to transversal current fluctuations. In section |^ we formulate a restricted not self consistent theory, which still 
contains all hydrodynamic modes and the most important rotational degrees of freedom. As a further, very drastic 
approximation, we set up a simple Maxwell theory to demonstrate that all qualitative features, observed in light 
scattering experiments in supercooled liquids are reproduced already within the restricted theory, by changing only 
the time scale of the structural relaxation. In section ^ the relation to other light scattering theories is discussed. 
The phenomenological equations of [l^ are derived within our theory and shown to be a special case of the restricted 
theory. 

II. LIGHT SCATTERING FROM MOLECULAR SYSTEMS 

Following p^ ] the total polarizability of a molecular liquid can be expressed as a sum of single particle and interaction 
induced many particle contributions. 

cx{r - r')E''''\r') = ^a'5(r - r') 5[r' - r')E'"'{r') 

i 

+ ^a'J'(5(r - r') Sir' - r^)E''^\r') 

ijk 

Where the superscript ext denotes the external electric field and int denotes an effective internal field. In this work we 
are concentrating on the direct contribution for linear molecules only. In a dense liquid other scattering mechanisms 
like direct and indirect dipole induced dipole scattering mechanism and collision induced scattering mechanisms are 
also present, in general. For CS2 e.g. it has been argued by Madden and Tildesley |^ that the interaction induced 
mechanism is the dominant one (although this result has been questioned For Salol it was found p8| , p^ , that 

the depolarized light scattering is dominated by scattering at orientational fluctuations. In more examples are 

discussed. A theory for depolarized DID spectra in simple liquids was presented in . A generalization of this work 
to a hydrodynamic theory of of light scattering, which incorporates in principle all possible scattering mechanisms 
was developed in [Q. Here we want to concentrate on the molecular origin of the depolarized light scattering for 
a system of linear molecules. A generalization to non linear molecules is straight forward, but to avoid unnecessary 
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complexity, we want to restrict ourselves to linear molecules. The numerical study of water |2U3^ has shown, that 
even the dynamics of this non linear molecule can be reasonably well reproduced by modeling it as as a linear molecule. 
Only direct contribution of the orientations to the spectra are discussed. All interaction induced contributions are 
neglected. 

For linear molecules the polarization tensor of a single molecule can be written in the form 

a« =al+|(37i(')n«-l) (2) 

Here h^'^'^ is a unit vector in the direction of the principal axis of the molecule i. I is the identity matrix. The numbers 
a and g are the mean polarizability and the anisotropy of the polarization respectively. 

The collective wave- vector dependent fluctuations of the total polarizability tensor a{q,t) are defined by 

N 

Sasi{q,t) ^^Sa^g]{t)e^p{iqri(t)). (3) 

1=1 

Here the indices S, I denote the direction of the scattered and incident light, respectively. The fluctuation of the 
single molecule polarizability dag\(t) is completely due to rotations of the molecule. 
The light scattering intensity is then given by 

Isi{q,u;) ^ FT{{a*sjiq,t)asiiq,0))) 



OO 



die 

— oo 



where the () brackets denote a thermal expectation value. Usually the convention is used that the xz-plane is the 
scattering plane and the scattering vector is anti-parallel to the z-axis (see Fig. |l|). In this geometry the depolarized 
scattering in VV and VH geometry, respectively, read : 

/^^(q, iu) = IjFT{{Sayy{q, tySayyiq, 0)))(c.) (5) 
/^^(q, uj) = IiFT {{5a;,{q, t)6ay^{q, 0)) sin2(e/2) 

+ {Sa*y,iq,t)Say,{q, 0))cos^{e/2)) {uj) (6) 

where G is the scattering angle, 9 = |<7| and FT denotes the Fourier transform. For linear molecules, the specified 
scattering frame is equivalent to the so called q -frame, where the z - axis is parallel to the scattering vector q, since the 
polarizability tensor is a spherical tensor of rank two, where only I = and I = 2 components appear (see appendix 
j). These are invariant under reflection q ^ —q. The explicit form of cx in the q - frame is derived in appendix |^. 



With (C6) we arrive at the result 



/^^(9,c^) ^ a25"°„ (9,^) + .92^(5^2(9,^) + y"UQ,^))-ag^^S"lo{q,u^) (7) 

/^«(<7,^)^g2^(sin2(e/2)5i2(g,a;)-fcos2(e/2)5"^2(g,c.)) (8) 
15 

In the derivation of Eqs. (0) and (^) we used that in the q -frame, the dynamic structure factors SJ^,™ {q,t) are 
diagonal with respect to m, m' and that for linear molecules S^,{q, t) = S'[ir^{q^ t). The contribution ^ *S"'2o('Z, w) can 
be neglected for light scattering experiments, since it is of relative order q^, compared to S"22iq,t) and S"QQ{q,t). 

Similar equations were already derived in |22| for the special case of a diluted gas, where translational and rotational 
motion can be factorized. In the case of dense liquids, the form Eqs. (|^) and (||) have to be used. For a system of 
linear molecules with only direct but anisotropic scattering mechanism, they provide a complete and exact description 
of the observed light scattering spectra. Since we are interested in light scattering experiments in supercooled liquids, 
we are going to present a microscopic theory for the generalized dynamic structure factors S";™ (g, t) which contains all 
observed features of light scattering experiments close to the glass transition. Already from the general form of (|^) 
and (^) some general conclusions can be drawn. 

First, due to the scattering angle Q in the VH geometry a backscattering geometry only observes the m — 2 
component of the dynamic density correlation function whereas a 90 degree scattering angle probes a mixture of 



3 



m = 1 and m = 2 components. In the following sections we show that the 5*22(9,'^) and 6*22 (f?,!^) component couple 
to the longitudinal and transverse sound mode, respectively, while the S22{q,^) component does not couple to any 
hydrodynamic mode. Therefore it is possible to replace 5'|2(g,a;) by its value at g = 0, but not the correlators 
with m — 0, 1. The different behavior of the correlators with different m is due to a dynamic breaking of the 
rotational invariance on the spatial scale of the light scattering experiments, caused by the existence of hydrodynamic 
singularities. The phenomenological form for the light scattering spectra used in the literature can be recovered by 
rewriting Eqs. and (^) in the form 

J^^(g,c.) ~ a'S"l„{q,cj)+g'^{P"Uq = 0,c.) + 1{S"Uq,^) - S"Uq = 0,u;))) (9) 

J^^(g,a;) ~ g'^iS^'^Q = 0,uj) + cos^e/2){S"Uq,cu) - S'l^ - 0,c.))) (10) 
15 

Here we have used that for q = due to rotational symmetry all correlation functions with different helicity m 
but the same / are equal and are therefore replaced with m = 0. In addition we have neglected the off diagonal 
contribution 6*20 (q, w). The last two terms in Eqs. (||) and ( p^ contain hydrodynamic poles and thus essential for 
the understanding of light scattering experiments. For theoretical considerations it is more convenient to use Eqs. (0) 
and (^). Therefore we continue in this work using the representation in spherical harmonics. 

Second, the spectra fulfill sum rules. The total intensities I^J^ — J du!lfj^{u!) for g — > are derived from Eqs. (|^) 
and (|) 

irj{q^O,t = 0)^a'S"°o{Q = 0,t = 0)+g'^ ^5"^2(9 = 0,t^0) (11) 

iriiq^O,t^O)^g^^S"Uq = 0,t = 0) . (12) 

Here we used, that the static structure factor Slj,{q = 0,t — 0) are independent of m and diagonal in I, I' . Therefore 
the anisotropic part of lYJ[ is exactly equal to ^lYo? ■ 

Third, we note the well established fact, that the often used procedure to obtain the isotropic scattering contribution 
(e.g. in |^3|, see also [^) by simply subtracting |/^^ from I^^ {uj) is in general not exact, since there are big 
qualitative differences between the correlation functions S^{q,uj) for different m, due to the coupling of translational 
and rotational motion. In Appendix 7B of |2^] the mentioned relation between I^^ {uj) and I^^ {to) could be derived 
by explicitly assuming, that rotational and translational motion are independent. But it is clear, that this cannot 
hold in dense liquids, where the rotation of a molecule can easily cause the build up of local stress via interaction 
with its neighbors. This stress can then be released by a center of mass motion of neighboring molecules. Instead of 
obtaining only the isotropic contribution the mentioned subtraction method will yield the following expression. 



,2 c"0 



a 



S"ooi'l,^)+ (13) 
<7'§(^"22('?,^) - ^"22(9 = 0,^) - 4cos2(e/2)(5"^2(9,^) " ^"22(9 = 0,..))) 

Eq. (^3|) shows, that even in back scattering geometry deviations from the purely isotropic scattering are to be 
expected. To demonstrate this point in more detail, we plot in figure |^ the value of the depolarization ratio in 
backscattering geometry (6 = tt) for a simple model which we derive from our equations (see section |^). To further 
demonstrate the result of Eq. (13) we have (using the same model discussed in section ^v|) further plotted the quantity 
5"22(q, w) — S" 22{q — 0, ti^))/ S"'qo('?, lo). This is done on a linear frequency scale in Fig. |l^ and on a logarithmic scale 
in Fig. Qq. It can be clearly seen that deep in the liquid phase (r = 1) it is practically zero whereas strong deviations 
occur esspecially around the Brillouin lines when supercooling the liquid. From Eqs. (^) and (^^ it follows that the 
depolarization ratio in general is given by 

S"l2{q,iu) + \S"l2{q,^) 



sin\Q/2)S"l2{q,u;)+cos^{Q/2)S"\2{q.Lo) 



^4^1 S"l2{q,^)-S"l2{q^Q,^) 

3 3y'^2(9 = 0,c.) + cos2(e/2)(5"i2(9,^)-^"22('7 = 0,c.)) 

In backscattering geometry (0 = tt) it is nearly constant in the liquid and given by 4/3 whereas strong deviations 
occur around the Brillouin line when super-cooling the liquid. These deviations are expected to be even bigger when 
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the experiment is performed not in backscattering geometry since in this case the m = 1 component which couples 
to the transverse phonon will also contribute. Fig. y shows the depolarization ratio on a logarithmic frequency scale. 
The deviations of the depolarization ratio from 4/3 are also seen in the experiment and are discussed in Therefore 
we have proven that the quantity lyv ~ 4/3/^// is not proportional to Sqq {q,uj), especially not around the Brillouin 
lines. We are aware that this is a severe difficulty in interpreting light-scattering experiments and want to point out 
that in the low frequency regime it is usually not applied anyways since lyn becomes very weak. 



III. THE EQUATION OF MOTION OF LINEAR MOLECULES 

After having an exact expression for the scattered intensities (see Eqs. (0) and (||)), we turn in this section to 
the second and independent step in order to arrive at a microscopic understanding of light scattering spectra. This 
second step is the discussion of a microscopic theory of a molecular liquid. The mode coupling theory of the glass 
transition has by now proven to be a very successful microscopic theory to describe the dynamics of supercooled 



liquids close to the glass transition (for a recent review see e.g. |24 25|]). In its early version |£6| only simple liquids 
consisting of spherical molecules were described microscopically. Due to its successful application also to experiments 
and simulations at arbitrary molecules it became necessary to develop generalizations to molecular liquids ■ 
Also aspects of these generalized theories were successfully tested in simulations [^0|-^. The main quantity, which 
is studied in is the coherent dynamic structure factor S]^,{q^t) for linear molecules. The theory describes on a 
microscopic basis the coupling of translational and orientational degrees of freedom, but does not take into account 
the coupling to transversal currents. In pT[ | the theory for the dynamic structure factor »S'Jp™„,((f, i) of arbitrary rigid 
molecules was derived. This theory deals also with the transversal currents, but is at present difficult to reduce to 
linear molecules. Since all theories of light scattering agree in the importance of transversal current fluctuations, 
it will be necessary to formulate a theory for linear molecules, which also contains coupling to transversal currents. 
Fortunately this theory is a straight forward generalization of the theory in Q . We do not take into account energy 
fluctuations in our derivation, which were necessary to describe the Rayleigh peak i.e. the influence of heat diffusion 
on the light scattering spectrum. We do not expect important changes for the discussion of the light scattering 
spectrum, by neglectingthc influence of heat fluctuations. For a discussion of the interplay of structural relaxation 
and heat diffusion see [ p3|jl4|] . 

The MMCT is derived within a Mori Zwanzig projection operator formalism for tensorial densities pim{q,t) 

N 

pim{q. t) = y/A^i' Yi^mt)) e^«~'(*) (15) 

with Yi„i{^) being the spherical harmonics and tensorial currents {j^yimiq, t). The index a — T,R denote translational 
currents and rotational currents, respectively. For linear molecules it is not necessary to use the orientational current 
as a vectorial quantity. The only quantity which will appear in the theory is 

rj- N 

{j^wm = ^/ Y.(^mmmt)) e^^-^^-w. (le) 

Here L is the angular momentum operator and uj is the angular velocity. {j^}im {q, t) with m E {—I, — / + 1, . . . ^ — 1, Z} 
are the components of an irreducible spherical tensor of rank I. 

To be able to describe light scattering spectra we will need all components of the translational current fluctuations 
{jji}im{q,t), n G {—1,0, 1}, not only its projection along q as in Here we used spherical components p, instead of 
Cartesian components {x, y, z}. 

N 

{jl}Uq,t) = z'V4^5]^;^ll™(a(0) e*«-^~'W. (17) 



i=l 



with 



{jli)lm{q.t) = -l{±{jl}i^{q^t) + l{j'^}i^{q,t)) 

blhmiq. t) - {]jAhn{q, t) (18) 
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The translational currents {jj^}im{q,t) can be written as a sum of components of an irreducible spherical tensor 

{J^hmiq, t) = J2 Caih ^lmm){f}Y;;iq, t) (19) 
l,rh 

with \l — 1\ < I < I + 1 and rri = m + fi, due to the properties of the Clebsch Gordan coefficients C{lil2l; mim2m). As 
an important special case we note, that the transversal center of mass currents {j'^i)Qo{q,t) have irreducible spherical 

components {j^}i'±i(9, i)- 

The currents and the densities fulfill a continuity equation 

^^Plmiq, t)^J2 'i-^r<l-,(l){j^}lrn{q, t) (20) 



dt 



vifri) m 



with 



With the help of Mori - Zwanzig projection operator techniques [[34yi[| , it is possible to derive a formally ex- 
act set of equations for the correlation functions S'j™™ {q, t) of the tensorial densities p;,„ (q, t) . For simplicity we 
write the equation in the q - frame i.e. {q = 9(0,0,1)). For a short ranged potential every correlation function 
{Aim{q,t)\Bi,„i>{q)) = {Aimiq,t)* Bvm'{q)) of two components Aimiq.t), Bi^rn'{q) of irreducible spherical tensors 
fulfill the relations 

{Airniq,t)\Bi,m'iq)) ^ S„,„,,{Ai.miq,t)\Bi,rn{q)) (22) 

{Airn{q,t)\Bi,rniq)) ^ ql'"''' for q ^ 0. 



The first line or (22) is just due to the fact that under a rotation around the z - axis, {Ai„i{q,t)\Bii„ii{q)) will 

transform into e*'™ ^'^'>'^{Aimiq,t)\Biim'{q))- Since the correlation function has to be invariant under this operation, 

the relation (^2|) follows. The second line is a consequence of global rotational invariance of the system of linear 

molecules. In addition, the correlation function {Aim{q,t)\Biimiq)) is independent of m for q = 0, if there are no long 

range correlations, which destroy the global rotational invariance. For correlation functions which behave regular at 

q — i.e. which do not contain any hydrodynamic poles, we are allowed to replace them for small q by their value 

at g — plus corrections. If there are hydrodynamic poles, the differences between different m can be crucial e.g. 

we will see, that the dynamic correlation function S^{q, z) contains for m = 1 even at small values of q couplings to 

transversal current fluctuations, where the one for m = 2 behaves regular at g — > 0. 

The dynamic correlation function S*;™™ (9,0 is a real quantity Q. Since we will discuss in the following part of 

the paper also current-current correlation functions we use the notation S™,{q,t) = (0 )J^™(g, t) The equation for 

=pp 

(j) {q, t) is therefore 
=pp 

d 

—6 {q,t)^-i q<j)^ {q,t) = -i (a, 

Ot=pp ==3P ^-^=0 =3Sp 

a 

^J>. {q.t)^~iT. {q)<t> {q,t) 
Ot=Jp —IP =pp 



I dt' M,M~t') <t> {q,t) 

Jo " =3 P 



^6 Aq,t)=iq^^^{q,t) = i yZq^cj) Jq,t) 

Ot=P3' ==33' ^=0 =3^3' 

^6 {q,t) = -iT. {q)<t> ^^{q,t) 
01=33' =P3 

- I dt' M..„{q,t-t') cl) {q,t) 

Jo — =3 y 



(23) 
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where we use the short-hand notation j = {j"}zm, j' = {jZ}i'm' and q = Siii{—l)'^q"_ (I), cj) il^i) = 

<5mm'+M({iw}i™('?;i)|Pim'('Z,0)), (q, = 0+ ((?, and = ^p+m.^'+m' ({in limC?, OI{iu''}i'm' 0))) 

are the current -density, density-current and current-current correlation functions, respectively. The matrix Fj (q) = 
<5/j,o<?"'^^(S^^ (g, rn)) is determined by the static molecular correlators. Qt and O/? are the mass and inertia, respec- 
tively. The products q^T_ .j^ (q) and g^F ^ (q) are the matrix of rotational and translational microscopic frequencies, 

=0 —lo P =0 —Jo P 

respectively. The memory matrix Miq^t) = M^^,' = {M;;;'}^';^,(g; = {QC{{]'^{q)}\R' {t)\QC{{]'^, {q)})^ is 
a frequency dependent damping matrix, C is the Liouville operator, Q is the projection operator perpendicular to 
the density- and current fluctuations (see Eqs. ([T5|), ( p^ and (|T7|)). R'{t) is a reduced time translation operator 
|. The memory matrix is not diagonal in ra and to', contrary to the one appearing in 0. This 
is due to the fact, that the currents for /i 7^ are not components of an irreducible spherical tensor. Instead, due to 
( |l9| ) and (|2^) the relation m' — m + 11 has to be fulfilled: 

Without the memory matrix equation, ( ]23| ) would describe a system of coupled undamped harmonic " modes" , where 
the modes are in this case correlation functions of tensorial densities. The physical origin of the memory matrix is the 
damping of these oscillatory modes due to interaction between them including translation rotation couplings, caused 
by the anharmonicities of the microscopic interaction potentials. Of special importance will be the induction of a stress 
^{jT}oo{q,t)/q by the force C{jQ}2i{q,t) caused by the rotation of the molecules. This mechanism is responsible for 
the existence of hydrodynamic singularities in auto correlation functions of non hydrodynamic fluctuations. 



R'{t) = Qe'Q^Q* 1 34 



A. Molecular mode coupling theory 



Within the MMCT, the memory functions are written as a sum of bare Markovian damping terms plus mode coupling 
terms. The mode coupling terms have the form of self consistent statically renormalized one loop approximations. 

K,^' (9; ^ - = » K",'' -t') + ^K^' (9; t-t') (25) 

The derivation of the mode coupling approximations is analogous to the one in Q . For the memory functions with 
fi — n' = the final result is identical to For general it can be written 

{<;:}rr'(.;o-^(£)'E'EEE>< 

qiq2 mim2 I1I2 I'^l'^ 

><{^^'}Di;"S('?,gi,92)5[:^^(gi,i)5S^(g2,i) , (26) 



with 



EK}r4t''"'('?'9i'92)CJ<7i) + (-1)"(1 ^ 2) 



(28) 



where cY\fiq) is the direct correlation function and 



r CKT m,ml,m2/ \ 



(2^i + l)(2/2 + l) 
{21 + 1) 



[l + (-l)'i+'^+'l 
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X 



qi{p)C{hl2l;000) ; a = T 



y/h{h + l)C{hl2l;101) ; a = R 
Here the functions qi{fi) are given by 



(29) 



At: 

91 (a*) = ^lyy (/^^^iMie^i.^.J + (1 - ImI) ilo(e,,,$,j) (30) 

The functions dl^,^{e) are related to Wigner's rotation matrices (we follow the notation of Gray and Gubbins 
p5|). For given Euler angles $, 8, x they are defined as |Q 

7?^„„,($,e,x) = e'""*4„,,(e)e-*™'>^. (31) 

,"I>q. are the standard spherical coordinates of ql with respect to q. The prime at the first summation in 
Eq.(p6|) restricts 91,(72 such, that qi + q2 — q in order to fulfill momentum conservation. Eqs. (p3|), (|26|), ( |27| ) form 
a set of self consistent equations for the generalized dynamic structure factors S]^,{q,t) of linear molecules. They 
are slightly more general than the equations in [Q| by including the coupling to transverse current fluctuations via a 
rotation - translation coupling. It is to be expected that this coupling will affect the results for the glass transition 
temperatures and the non crgodicity parameters studied in [l|j30[| only quantitatively but not qualitatively p6| . The 
dynamics, instead, can be changed qualitatively in certain wave vector - ranges. Especially for small wave vectors the 
hydrodynamic pole in the transverse current fluctuations can have large effects on the density relaxation spectrum. 
We will demonstrate explicitly in the next chapter, that the coupling to transverse current fluctuations is necessary 
to reproduce the appearance of transverse sound modes in Brillouin scattering spectrum of linear molecules, within 
the framework of MMCT. 

To calculate the light scattering spectra it is most convenient to perform a Laplace Transform of Eq. ( p3| ) . With 
LT{f[t)){z) — i Jq°° e*^*/(t), with Im{z) > 0, we obtain the following matrix equation 




(z) M 



=PP =P0 — _ =PP =P3 

0. [z) (/)..(z) ) " I ^0 ^0 
=jp =jj / \ =jp =]j 



(32) 



Here we have chosen a simplified notation. In the q - frame the first matrix in Eq. (^2|) would be explicitly: 

/ 4 -g zSu'S^„,, -<5n'5™™'(-l)"g%(0 A 

V -Ejp 4 + MjJ ; V -^^,0 ^qS 5^m' (S-i)((7) z Si, 5^„, + M^^;, (g, z) ) ^^""^ 

The matrix of static correlators on the right hand side of Eq. (5^) is 





















\ =3P 


=J3 / 





(34) 



It is obvious from equations (32) and (^3|) that the sum rules ( |lT|) and ( [l^ ) are automatically fulfilled, if the 

approximations for the memoryfunctions obey the very weak requirement, that lim^^^oo+ie ^ hm^^"' — _ where e is 
an arbitrary positive number. Since memory functions are regular for t = in most physical cases, they even vanish 
as l/z for z ^ 00. This property is especially fulfilled for the selfconsistent MMCT and also for the primitive theory 



used in section (iV) 



IV. GLASSY DYNAMICS AND HYDRODYNAMIC MODES 



Light scattering usually measures at small wave vectors q. Therefore the correct treatment of hydrodynamic modes 
becomes crucial for explaining light scattering experiments. To obtain the generalized density correlation S^,{q,t) for 
small wave- vectors, it would be necessary to solve the self consistent set of equations ( |2^ , |26| , p7| ) for all wave- vectors and 
all l,m, since all degrees of freedom are coupled via mode couphng integrals. To study the hght scattering problem, 
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we first want to restrict the discussion to the fluctuations, which are most relevant for the understanding of hght 
scattering experiments. These are the density fluctuation pim{<l, t) for (Z, m) = (0, 0) and 2, m their respective current 
fluctuations {jo}ooi<l,t), {jo'}2rn{q,t), which are directly measured in light scattering experiments (i.e. I = 0,2 and 
m = 0,1,2) and the current fluctuations {iJ}oo('?,0; which are slow, since the total currents {jJ}oo(<7 — 0,t) are 
conserved. Here we have also neglected the current fluctuations {jo}2mi<l,t). Their contribution will be of higher 
order in q, as can be seen by comparing the translational and rotational current contribution in the first line of ( p^ ) 
together with (|2^). Additional simplification occur in the limit q ^ 0, due to relation (|2^) for correlation functions 
of spherical tensors in a rotationally invariant system. 

1. The static correlation function is Sll]{q = 0) = Sni Si = ^["(g = 0), independent of m. 



2. The component of the matrix F,^ of Eq. ( |32| ) with I = I' = m — fi — and a — T reduces to — c^\<l- 
Here cy is the longitudinal isothermal sound velocity in the liquid. 

3. The component of the matrix Fj^ of Eq. (|3^) with I = V = 2, ji = Q and a ~ R reduces to e^(g) ^ ~ 
where cufj is a classical frequency, related to the rotation of the quadrupoles. 

4. The memory matrix {M(^"^}gQ(iO) has to reduce to irjiq^ , where rji is the longitudinal viscosity of the liquid. The 
transversal memory function is for q — > 0, z — > «0 given by {M^^}QQ{iO) = iq^rjs, where 775 is the shear viscosity. 
The q^ dependence of the translational memory functions {w^J/}og are due to momentum conservation. The 
parameter {i^JJ^jgo and {tJJ^Iqq in are identical to the longitudinal modulus Ki and the longitudinal a 
- relaxation time t; for /i = /i' = and the shear modulus Gs and the transverse a - relaxation time Tg for 
fi = ^' = 1, respectively. With this choice the Maxwell relations iji — 7]° + Kitb and 775 = 77^ + Ggrg, where 
Tjf, rig are the contribution from the Markovian part of the memory matrix, arc fulfilled in the liquid. 

Due to the local nature of the cage effect, which is responsible for the slowing down of structural relaxations, it is 
strictly speaking not possible to study self-consistently the hydrodynamic limit, without the knowledge of relaxations 
on local length scales r oc rg, where rg is on the scale of intermolecular distances. But since the memory functions 
m"'", (q, t) do not contain any hydrodynamic pole by construction, they are non trivial only due to the glass transition 
dynamics, which in turn is independent of the hydrodynamic fiuctuations at short wavelength. It is e.g. theoretically 



understood |37|] and verified in simulations |38|, that systems with qualitatively different hydrodynamic behavior 
exhibit the same glassy dynamics. To obtain the qualitative behavior of Eqs. ( p3| , p6p7| ) , it is therefore sufficient to 
replace the mode coupling part of the memory function matrix by its leading wave-vector behavior multiplied with 
a function, which is able to describe glassy dynamics. Although the memory functions are free of hydrodynamic 
singularities, they exhibit the full frequency dependence of glassy dynamics ("fast" and "slow" (3 - relaxations, a - 
relaxation plus additional complications as e.g. contributions from Bose peak phenomena [ ^9pO| ]). For reproducing all 
details of light scattering spectra, which are directly related to glassy dynamics, two approaches are possible. Either 
an ansatz has to be foun d, w hich is compatible with all the mentioned phenomena, or the full set of microscopic 



equations, derived in sec. Ill, had to be solved numerically, to account at least for /3 - and a - relaxation, (see e.g. 



yj). But before this very difficult problem can be treated, it is necessary to demonstrate, that the structure of 



the equations derived in sec III can account for all the hydrodynamic poles and their interplay with the most basic 
phenomena of structural relaxations i.e the a - relaxation. To achieve that, it is sufficient to use an a - relaxation 
ansatz for the memory function. To make the analysis as simple and explicit as possible, we choose simple exponentials 
(Maxwell theory) for the non vanishing memory functions, multiplied with their leading wave vector dependence, q", 
n €{0,1,2}. 



The small wave vector behavior of the memory functions (|35| ) can be derived from ( p2[) and (|1S| ) and the conservation 
laws for total momentum in every spatial direction for I — or I' — 0. The relaxation times {t^'^> }™'™ are taken at 
q = 0. The projection operator formalism guarantees, that this value is nonzero, since the memory functions do not 
contain hydrodynamic poles. The Markovian part of the memory functions is in the following neglected, if it would 
vanish at wave vector q — 0. There are two severe consequences of this approximation together with the a - relaxation 
ansatzes for the memory functions. First the sound poles in the glass do not show any damping, instead of a damping 
proportional to q^ . Second, a fit with the Maxwell ansatz or any other ansatz, which only describes the a - relaxation, 
would lead to an artificial time scale separation for the a - relaxation times of the different memory functions (see 
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e.g. the discussion in The source of both errors, is the neglect of "fast" /3 - relaxation phenomena, which also 

contribute to the memory functions. 

The task is now to show, that the predicted spectra are consistent with sum rules and the qualitative behavior of 
light scattering experiments, which already follow very generally from a purely generalized hydrodynamic analysis 
p4| and that the interplay of rotational and translational motions lead to the qualitatively correct renormalizations 
of the hydrodynamic poles, when entering the glassy regime. 



If we order our basic variables in the form poo, P21 
in Eq. (^l are 



1, {jo }2in, {ii }oo, the frequency independent matrices appearing 



V6 



V6 








/41 



\/6 0/ 



\ 



/ 



(36) 



The slow part of the memory matrix is given by 



zl- 



qKiR Tin 



2 tirM 

Kr Tr 

Z TR + i 



\ 



. _ Kr Tr 
ZTR + i 

qKsR TSR 
ZTSR + i 



Kr, Tr 
z TR+i 



qKsR TSR 
z TSR+i 

_ q^GsTs 
ZTs+i 



(37) 



The leading wave-vector dependence for q — > of the matrix elements is derived from conservation laws and Eq. 
(p^). All the matrix elements, which arc left empty are exactly zero due to m^m! selection rules (see Eq. (p^)). For 
the appearance of transversal sound modes in the light scattering spectra and the explanation of the Rytov Dip it 
is crucial that Ksb. 7^ 0, which quantifies the memory matrix-element between the transversal current for I — and 
the rotational current for m — 1. The four parameters Ki, Kn, Km and have to be such that the memory 
matrix remains positive definite for all frequencies. For z = the relation 



KlKjl TlTR 



K'm T^R > and KrGs tstr 



KIr 



SR 



> 



follow and the diagonal elements Ki,Kr and Gs have to be positive. We also note that the exact relations 

GsKr > Kg]^ 

KiKr > Kfj, 

can be derived from the Cauchy relations 

by considering the low and high frequency limits of Eqs. ( p| ) and (^). 



(38) 



(39) 
(40) 



(41) 
(42) 



A. Hydrodynamic poles 

Before we discuss the numerical solution of Eq. (^^, it is useful to investigate its hydrodynamic poles. The only 
conserved quantities are (besides the total energy) the center of mass density and the total momentum in every 
spatial direction. They are the cause of the hydrodynamic poles (i.e. poles which show dispersion laws z oc g", with 
n = 1,2) in the respective auto correlation functions (<i>pp)oo and {($jj) JJloo}- But due to the dynamic couphng 
of the rotational degrees of freedom and the translational degrees of freedom, which appear naturally in the memory 
matrix, also the correlation functions of non hydrodynamic variables do exhibit hydrodynamic poles. To study this 
phenomena, we invert the matrix in Eq. (|3^), use the ansatz z = pq" , n e {0, 1, 2} and expand the denominator in 
powers of q. For n = this gives poles of non-hydrodynamic nature ( rotational modes with a frequency p), for n = 1 
propagating modes (transverse or longitudinal phonon modes) are described with p being the sound velocity and for 
n = 2 a diffusive mode with a transport coefficient (in our case generalized viscosity) —p/i is obtained. 
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Let's first study the transversal current fluctuations. In a simple liquid the transversal current correlator 
i(^jj)'ii"}ooi^) exhibits a viscous pole at z = —iGsTsq^. This is also the case for the liquids of linear molecules, 
studied in this paper. But in addition, due to the translation rotation coupling, also the auto correlation functions 
of densities pim.{q,t) or currents {i(f}im with I = 2, m = 1 do exhibit the transversal hydrodynamic poles. The 
reason for that is, that the tensor {jf}oo has irreducible spherical components with / = 2,m — 1, and thus is able to 
couple dynamically to the specified tensors via the memory functions with m = 1. There is also a non- hydrodynamic 
singularity related to the rotational motion of the molecules. If we restrict the correlators to their poles (neglecting 
the glassy dynamics), the poles of the ($pp)22(2), {(%j)ocf }22(2)) {(%j)i7}oo(^) components are given by: 

(z' + i{Krtr + yR)z ~ Lol) {z + iGs rgq^) = (43) 

The two poles couple into the dynamic correlators with different amplitudes. The coupling of the second pole into 
the I — 2, m — 1 rotational component causes the Rytov dip. As an example we therefore give the the term in lowest 
order of q of the strength of the transverse sound mode coupling into the ('I'pp)22('^) component. In the vicinity of 
the Rytov dip this component can be expressed as: 

Therefore the strength of the Rytov dip is proportional to q^ times the matrix element which couples the transverse 
current to the m — 1 rotational motion and vanishes if the rotational frequency ujr goes to infinity. 

As soon as zts ^ 1, the diffusive pole will turn into a propagating transversal sound mode. In simple liq- 
uids this pole will be at z = zLy/Gsq- Whereas in molecular liquids the transversal sound velocity is renor- 
malized by contributions of the rotational degrees of freedom. The pole structure of the specified correlators 
(*pp)22(^)> {(*Jj)Sf }22(^). {(%)?7}oo(^) is in the supercooled regime given by: 



(z' + ^VRZ - Lol - Kr) ( z^ - {Gs - TT^^^M' ) = (45) 



{Kr + u;1)' 



I.e. the transversal sound pole is given by 



MlGs-j^^^q :=±c,q (46) 



The transverse sound velocity is shifted to smaller frequencies compared to what is expected in a simple liquid. This 
trend was already noted in IQ. Note, that due to the positivity of and the exact relation (^, the transversal sound 
velocity is always well defined. I.e. by treating the rotation translation coupling explicitly, we are able to describe the 
contribution of the rotational motion to the transversal sound velocity c±. The microscopic rotational translational 
coupling is the cause of the appearance of hydrodynamic poles in correlation functions of non hydrodynamic (i.e. for 
q non conserved) variables and of a renormalization of the transversal sound velocity. Analogous behavior is 
found for center of mass and longitudinal current fluctuations. The amplitude of the transverse sound pole in lowest 
order of q which occurs in the (*I'pp)22(-^) component can be derived. It is in the vicinity of the transverse phonon 
mode given by: 

)i ( KliGsiKR + u^D) \ -1 .47^ 

{'^PPh2K-^} y^KR + u;l){KsR^{Lol-KR) + GsKR{KR + ul))) ^^'^ 



z±JG. 



^SR 



(Kr+uj-J 

Therefore the transverse phonon mode can only be observed as long as the rotation couples via Kr to the structural 
relaxation. 

In the fluid (TiCyq ^ 1) the longitudinal components with I — Q,m — Q and the rotational components with 
Z = 2, TO = contain two types of modes . First the I = 2 rotational mode and second the longitudinal phonon mode. 
The poles of (<i>pp)Oo(^), {^pp)U^). {(%)?o^}§^(z) are given by: 

{z^ + i{KrTr + vr)z ~ ujI) _ c|(?2 + iRmzq^'^ = (48) 

The first pole gives the damped rotation of the molecule, whereas the second term describes the usual longitudinal 
sound modes in the liquid at z± — icng with the damping Kinq^ /2. This expression is vaHd as long as Kinq^ <C 
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TiCi\q ^ 1. The amplitude of the longitudinal phonon mode in {^pp)22{z) is proportional to g^, where it is of order q° 
in pp)qo{z) ■ Its contribution to the sound pole in the - spectrum in the liquid can therefore be neglected. 

In the solid {Tic\\q 3> 1), the sound pole will be shifted to higher frequencies and, as an artifact of the Maxwell 
theory, the damping vanishes. An inclusion of (3 relaxation phenomena will cure this unphysical behavior. We obtain 
as poles of the ($pp)go(^), (*pp)^2(^), {(%)Jo^}oo(^)> {i'i'jj)m }22{^) components as: 

(z^ + zi^nz -ul- Kr) (z^ - (cjj +Ki- )qA - (49) 

i.e the longitudinal sound velocity Coo is, as the transversal sound velocity, modified by rotational degrees of freedom. 

cl^cj + Ki- , (50) 

Due to the positivity of ajj^ and the exact relation (|4^) the sound velocity is always shifted to higher values in the 
glass, but the shift is reduced compared to what would be expected in a simple liquid. In analogy to the transverse 
mode we can give the low q expansion for the amplitude of the longitudinal sound pole in the (*I'pp)22(^)- the 
vicinity of the longitudinal phonon frequency this component is given by: 

M22{^) = (iF^) ; (51) 

zTyC|| +Ki- ) 

Therefore a longitudinal phonon in the (*I'pp)22(^) component is always observable in the supercooled liquid as long 
as the rotational motion couples via the matrix element Kr to the structural relaxation. 



B. A solution of the equation of motion 

We have now solved the equation of motion for some chosen but fixed parameters. Close to the glass transition it is 
only the scale of the a - relaxation time, which is changing considerably. If the time temperature superposition principle 
(TTS) is fulfilled exactly all relaxation times {r;^^/}^™ (T) can be written as {r;jp'}jp™7(?')- The function /(T) is 
the same for all relaxation times and is a quickly changing function of temperature. In the idealized mode coupling 
theory it would be {\T -Td/T^)-^ , often used fit formulas are the Vogel Fulcher function /(T) = Aexp{-B/{T-To)) 
or as a special case the Arrhenius law {Tq = 0). The prefactors {t^^/ are constant, but in general different for 

different combinations of jj! ,a,a' ,1,1' ,m,m! . Very often the TTS is violated at lower temperatures. In this case 
also the {f^^, }^™ vary slowly with temperature. 

Since we are in this paper only interested in qualitative aspects of the solution, all a - relaxation times were taken 
equal i.e. {r^^/ }IF'™ = 1- The function f{T) is called r in the following. In this way we demonstrate how we 
can explain the complete viscosity range from the highly viscous liquid towards the glass by just varying a single 
parameter, r. 

In the figures ^- |l^ we have chosen the following parameters: The frequency scale for the rotational frequency was 
set to unity = 1. In these units the other parameters where chosen as cy = 0.6, Gs = vr — 5, Ki = 1, Kr = 1, 
KiR = Ksr = 1/2. The external momentum q selected by the scattering experiment is set to q — 0.02. For zt <C 1 
the scattering experiment probes an amorphous solid whereas for zr ^ 1 it probes a liquid. 

In Fig. ^ we have plotted from a solution of Eq. ( p2|) the imaginary part of the center of mass correlator 5*00 (<?, z) 
which according to Eq. (^ forms the main part of the polarized light scattering intensity /^^(w). For r = 1 <C (cyg)"^ 
the Brillouin line, caused by a longitudinal sound wave, is at c\\q = 0.012. The damping is proportional to q^ (see 
Eq. (^)). For r = 100 ~ (cyg)"^ a broad central peak occurs together with a strong damping of the phonon modes. 
This is a well known phenomena which was already explained by Mountain for the Brillouin spectrum. The same 
mechanism applies here. If we calculate the hydrodynamic sound pole with the condition lot = 1 the equation for the 
sound pole is 

u;^ + icj^q^ + t^q^ - {cl + ^)q^ = (52) 

This equation can be solved with the ansatz w = Cooq + iTq + 0{q^) i.e. the damping is of order q ^ q^ instead of 
order g^. This eff'ect is even stronger in real glass- formers, where the a - relaxation is better described by a stretched 
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exponential behavior instead of a single exponential. We can estimate, that for frequencies in the high frequency wing 
of the a - relaxation, where the von Schweidler law applies for the memory-functions (m(z) oc (— iz)^''), the pole 
becomes a cut at = ±(coo + ir)q^^i . If we include (3 - relaxation like phenomena with Loni" {lj) ^ uj"", we obtain in 
the frequency range, where this fractal behavior holds, two strongly damped propagating modes at w = ±Coo9+*rg^"'"'^. 
Since the fractal behavior is experimentally observed in depolarized spectra even below Tg, where the a relaxation is 
far below the experimental frequency range (see e.g. [|3|), it could account for the anomalous strong damping of the 
Brillouin line in the sub Tg regime p3| . In this context it is important to note, that the actual physical reason for the 
appearance of a fractal part of the susceptibility spectrum is not relevant for the described mathematical mechanism 
of producing strong damping of phonon modes. The f3 - relaxation phenomenon close to the critical temperature of 
mode coupling theory is as good a candidate as the still not yet understood fractal behaviour below Tg ( p3|-^). 
For T — 10^ S> (cyg)"^ the light scattering probes a solid with a well defined, now renormalized, phonon mode with a 
renormalized sound velocity which is given by Eq. (pO[). As discussed above, the speed of sound in the glass is always 
bigger than the speed of sound of the liquid. 

The next figure^ shows the spectral function of 522(9,2), which gives the only contribution in forward scattering 
direction and zero contribution in backscattering geometry for the depolarized geometry due to the cos^(0/2) factor 
(see Eq.(||)). The 522(9,2;) component is strongly influenced by the coupling of transverse currents to the m = 1 
rotational current. 

iQC{Ai{q)\—^\QC{jli}ooiq)) « (53) 
z - QCQ ZT + I 

Deep in the liquid for r = 1 ^ {c\\q)^^ the purely diffusive transverse currents cause the Rytov dip at zero frequencies. 
This is shown in the inset of Fig. |^ When super-cooling the liquid the Rytov dip disappears and a broad central 
peak shows up for t = 100 w [c^)^^ which develops a shoulder at the frequency of shear waves. When further 
super-coolin g th e liquid for r = 10^ (cng)^^ the propagating transversal phonon modes of the solid at ±c^q shows 



up (see Eq. (|46D). 

A further contribution to the spectrum of depolarized light scattering according to Eq. (|^) is the m — 2 component. 
Due to the factor sin^(0/2) it is the only contribution in backscattering geometry. Since it does not couple to any 
hydrodynamic mode, it has the simple form 



S 22 (cc') = 5^ (9 = 0) , V r- (54) 



The spectral function of 512(9,2) in Fig. || shows for r = 1 <C (c^q)^^ no structure for low energies (far below the 
rotation spectra). It has no Rytov dip since there is no matrix clement which couples shear waves to the m = 2 
rotational currents. The absence of a Rytov dip is shown explicitly in the inset of Fig. ||. When super-cooling the 
liquid a broad central peak shows up for r = 100 w (c„9)~^ which narrows in the solid for r = 10^ ^ {c^q)~^ . 
In our units the height of the peak is always Kj^t -\- vji and the width is of order 1/r. The absence of transverse 
modes in backscattering geometry is clearly seen in experiments (see e.g. p7[ or earlier where light scattering 
in backscattering geometry was applied to the molecular glass former salol). In other scattering geometries the VH 
- spectrum exhibits ideally only the transverse phonon (for luts S> 1). We want to point out that there are two 
qualitatively different mechanism which lead to the observation of a phonon line in a depolarized light scattering 
geometry. The first one is leakage of the longitudinal phonon mode due to an imperfect polarization filter and should 
not be present in an ideal depolarized geometry. The second one at a lower frequency is a direct transverse phonon 
which couples to the 522(9, w) component. 

A further contribution which we have plotted is the 522(9,1-^) component. It enters according to Eq. (|^) into the 
intensity for the polarized light scattering geometry. There is a non zero matrix element which couples the longitudinal 
phonon mode to the m=0 rotational mode which we approximate from Eq. ( |3^ ) as 

{{jUo^{q)CQ{Q{z - C)Q)-'Q£{j,f}2,{q)) « (55) 

Therefore the longitudinal phonon mode shows up in the m — component when the liquid is supercooled. This 
is shown in Fig. For t = 1 no coupling to the longitudinal phonon mode can be observed. This can easily be 
understood. The coupling term ( p5| ) reduces to iKiETq for c\\qT <C 1 at the position of the phonon mode uo — C||9. 
Since the coupling term to the phonon mode will appear in second order perturbation theory in 9, the phonon mode 
is multiplied with a factor {KmTq)^ i.e. the maximum height of the phonon contribution in the spectrum 5"22 will 
be of order Kf^r/Ki, which is equal 1, in the units we are using in Fig. |l^. For Cj|9T » 1 the couphng term ( |55| ) 
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is Kiq/id i.e. at = c\\q, the coupling constant is of order 1 and the height of the phonon mode is again of order 
Kf^r/Ki ^ 1. The phonon mode can be detected as soon as c\\qT ^ 1. In our units this happens for r ~ 50. In Fig 
|lO| a broad shoulder can be seen for r = 100, which turns into a clearly defined phonon mode (r = lO'', r = 10^). 

The S"QQ{q,uj), S"22{q,^) and S"'22('Z, w) sum up according to Eq. (0) to give the total polarized hght-scattering 
intensity. This is plotted in Fig. |l^ where we have plotted a fictitious VV-spectrum under the additional assumption 
that the squared isotropic part of the polarizability is ten times as big as the anisotropic one (a^ — lOg^)- Note 
that the orientational correlator S^{q,uj) and the translational center of mass component iS'go((?,w) are of completely 
different origin even though a mixing of the poles occurs in the supercooled regime. This is best seen in the liquid 
regime where the microscopic frequency of the orientation is a rotational motion whereas the microscopic frequency 
of the center of mass component is given by the longitudinal phonon mode. It is experimentally impossible to extract 
S'qq((7, w) out of a measured spectrum except for the case of vanishing anisotropy. 



C. Susceptibility spectra 

In order to make the influence of structural relaxations more transparent, we have plotted on a logarithmic frequency 
scale the spectral functions of the susceptibilities 

x";r'(9,^)-^^";r'(9,^) (56) 

which correspond to the quantities plotted in Figs. |^, ^, ||, |l^ and ^ Fig. ^ shows the spectrum x"oo('?i^)- The 
central peak of Fig. ^ turns into an a-peak which upon super-cooling the glass transition moves out of the microscopic 
frequency which is given by the longitudinal phonon mode. We want to stress again that due to our particular simple 
choice of the memory function (Maxwell theory) the a-peak shown in Figs .|5|,|7|,p|, p]jT^ does not have the correct 
stretched form known for glassy systems and there is no true (3 relaxation regime. If the molecular system has no 
anisotropy (in Eq. (0) = 0) this contribution to light scattering is the only direct one. The area under the a-peak 
la on a logarithmic frequency scale, which is also the area under the Mountain peak on a linear frequency scale 
and the overall intensity It^t £^re in this case related to the non-ergodicity parameter (Edwards- Anderson parameter) 
f^i^{q,Lo) = la/Itot for the center of mass component. 

The figure (Fig. ^ shows the spectral function of X22('?j^)- The microscopic is given by the rotational I — 2 
mode roughly at wr. When the liquid is supercooled towards the glass transition the shear wave shows up when 
the a-relaxation moves over the frequency range range for transverse phonons at w = c±q. As expected the m = 2 
component x"22(.'lt'^) Fig. ^ only shows the broad rotational mode and the a-relaxation since there is no matrix 
element in the memory kernel which couples to the m = 2 component. 

Similar to the m = 1 component the m = component also shows a coupling of a hydrodynamic mode. However in 
this case it is the longitudinal phonon which couples into the m = susceptibility. The spectral function of X22('Z: ^) 
on a logarithmic frequency scale is plotted in Fig. O. 

The overall situation for a fictitious spectrum uW^ {q,uj) is shown in Fig. |l^ for the same parameters as in Fig. 
|T2| . Note that due to the mixture oil — V — and I — V — 2 components the Brillouin line of the longitudinal phonon 
consists of two components. One is caused by the direct observation of the center of mass component where the 
longitudinal phonon gives a low lying microscopic frequency whereas the other contribution comes from the phonon 
mode coupling into the I — I' = 2, m = Q component. 



D. Light-scattering near an orientational instability 

Further physics which is contained already in Eq. ( p^ are some aspects of light scattering near an orientational 
instability (e.g. near an isotropic-nematic transition). From Eq. (^) one can see that the amplitude of the light 
scattering intensity is for the depolarized light scattering spectra mainly given by the static density correlation S^{q). 
Close to a nematic transition, a weakly first order phase transition, S22{q) increases strongly for small q where 
limq^o '5'22('Z)~^ — HKerr is the optical Kerr constant. On the other hand the rotational frequency lu^^ (in Eq. (|3^)) 
contains S^{qy^^ in the denominator. In this way our equations describe the broad central peak together with the 
strong scattering intensities at the isotropic-nematic transition. 
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V. RELATION TO OTHER THEORIES 



The main issue of the paper so far, was to demonstrate that the theory of light scattering can be brought in a form, 
which is accessible to tested theories for the dynamics of supercooled molecular liquids and to show that the structure 
of the equations of motion reproduce the light scattering experiments for molecular liquids. We now will show how 
the phenomenological equations underlying the most recent theory of light scattering for molecular liquids can 
be rigorously derived within our theory. To keep the derivation as simple as possible, we will restrict ourselves to 
the variables poO: P2m, {jjlooj }2m used in chapter |^ We also will comment on and on [Q. In appendix ^ 
we demonstrate for the theory of Anderson and Pecorag how light scattering theories for linear molecules based on 
projection operator formalisms can be related to our theory. 

Dreyfus et. al start by writing the continuity equations for the center of mass density fluctuations and the 
center of mass momentum density fluctuations, 

d 

-Q^Poo{q,t) ^iqj{q,t) (57) 
d 

-g^ji{q,t) = iqj(Tji{q,t) , (58) 

where ji are the Cartesian components of the center of mass current fluctuations and aij are the Cartesian compo- 
nents of the stress-fluctuations. To obtain a closed set of equations it is necessary to write down constitutive equations 
for the stress-fluctuations which relate them to the current and density fluctuations. Instead of using phenomenolog- 
ical ansatzes, as it was mostly done in the existing literature, we will write down exact equations for the stress-tensor 
fluctuations using generalized constitutive equations. By using a formalism introduced in p3| we can express the 
stress-fluctuations for vanishing amplitude of the wave- vector exactly by the fluctuations of the basic set of variables, 
which we used in section |v| . (For simplicity it is more convenient to use the Cartesian components of the center of 
mass current fluctuations) 

cr^J = Poo (9,^) ^0^ J ^ (Poo|t,j) 

+ i]Mt)qk ® {(Jrk\R' {t)\a,j) (59) 

+ ^{jo}2ra{q,t) ® {C{j!^}2n,\R'{t)\a,,) 

Here, repeated indices are summed over and we defined = dt' A{t')B{t-t'). R'{t) is again the reduced time 

propagator acting in the space perpendicular to the chosen density and current fluctuations. Therefore no terms of 
the form (Cpim\R'(t)\'^ij) appear in equation (p9|). There are in principal terms proportional to the fluctuations p2m- 
But since in Eq. (p7h only the combination qjcrji appears, the proportionality factor contains the term qj{p2m\o'ji) — 
iP2m\^ji) = [C.p2m\ji) = i{jo}2m\ji) ~ . Therefore no fluctuations proportional to p2m do contribute to the 
generalized hydrodynamic equations. For g — > the equations ( p9[ ) considerably simplify. The only symmetric tensor of 
rank four {crrk\R'it)\(^ji), which does not vanish for g ^ is {arklR' it)\<7ji) = ip\R' {t)\p)6krSij + {a^^ \R' (t)\a'^j){SrjSki + 
SriSkj — "^SijSrk), where p = \ X], (Ja is the scalar part and afj the traceless part of the stress-tensor fluctuations. 
The third row of equations (B9h can be evaluated by transforming ct^ to spherical components aim- Then, for q — > 0, 

the identity {jQ}2m\R' iM^im') = <5i2^mm'('C{j^}2o|^'(i)k2o) and pim{t) = V^{j^}2mit) can be used. Therefore the 
constitutive equations ( p9| ) reduce to 

= (5F(t)(5y + ® iqj{t)5ij (60) 

+ -vs{t)<»n^{t) (61) 
n 

with ?7_B(t) = fceT V (-PI^'(OIp)' Vsii) — ^.^ ^ y (cr|j | -R' (t ) | crfj ) being the generalized bulk viscosity and shear viscosity, 
respectively ]3^ . The tensor — i{qijj + qjji — ^qj^ij) is the strain tensor. The fluctuations SP{t) of the internal 
hydrostatic pressure due to density fluctuations are given by KpQQ{q,t), where K = , is the (static) bulk 

module, i.e. the inverse of the compressibility n. The tensor Q is given by 
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and the function ii{t) is the the matrix element fi{t) — ^^^^ {C{jQ}2Q\R' {t)\T2o) . Now we only need another consti- 
tutive equation for the tensor Qij{t). Using the same strategy as in the derivation of the constitutive equations for 
the stress tensor we easily derive 

^Q.j [t) = -u;lQ,, {t) + ii{t) ® T,, (t) - {Mo«o«}°0(t) ® Q,, (t) (63) 

Here the same function as in Eq. ( pO| ) appears naturally within the formalisrn, confirming the Onsager principle. 
The memory- function {M^^}22{t) is the same as used in section [11. Eqs. (57) - ( |6^ are exactly the equations used 



|T|. With our formalism, we can identify the phenomenologically introduced functions fi{t) and Qij{t) of jl^. We 
also want to emphasize, that the convolution integrals in time i.e. the retardation effects are a necessary consequence 
of the slowing down of structural relaxations and its effect on the frequency dependent viscosities and the rotation 
- stress coupling function /i(i). The ansatz of Quentrec where the viscosities and /i(i) are replaced by only 
temperature dependent functions is therefore not acceptable for the description of supercooled liquids. 

It is important to note, that the form of the equations (|57| ) - ( p9| ) depends crucially on the chosen set of variables. 
If we would not have chosen the rotational currents {iQ}2m{<i^ t) explicitly as a member of our basic set of equations 
the last line of Eq. ( |60| ) would contain a coupling /t(t) ® Qij (t) to the tensor Qij instead of to its time derivative. The 
function jl{t) can also be expressed in terms of a (modified) reduced time propagator R{t), jl{t) cx ({j(f }2ol^'(OI'''2o)- In 
addition, the equation for Qij{t) were of first order in time instead of second order. Wang uses a mixed representation 
. His constitutive equation for the stress tensor coincides with Eq. (^) , but the equation for Qij is only of first 
order. For deriving exactly such a set of equations, it were necessary to use different projection operators for deriving 
the constitutive equation and the equation for the tensor Qij. From our point of view there are mainly two reasons 
why it seems more advantageous to choose one single set of basis variables including the currents {jo}2m{<l,t). First, 
since approximation schemes for force - force autocorrelation functions i.e. the memory-function {M^", (t) seem 
to be easier, than for mixed current - force memory- functions, which would appear, when only the densities pim{qit) 
and the conserved currents {jj}im('f, i) were used as variables. Second, by using one set of variables the Onsager 
relations are automatically fulfilled, since there appears the same function /x(t) in the equation for the stress and the 
tensor Qij{t). In the approach of Wang it is important to choose the approximations for the different functions jl{t) 
and fi{t) carefully, such that the Onsager principle is guaranteed. (Essentially the time derivative of ji{t) is related to 



The approach of 
to derive the genera 



L4| is more general than ours, since no assumptions on the scattering mechanisms were used 
form of light scattering spectra of supercooled liquids. This was achieved by only using the 
hydrodynamic variables center of mass density and center of mass currents as the basic set of variables for applying 
constitutive equations to the dielectric tensor fluctuations. In this way the coupling mechanisms between rotation 
and translations in molecular liquids are not explicitly treated but lead implicitly to frequency dependent Pockels 
constants, relating the hydrodynamic modes to the dielectric tensor fluctuations, and unknown background spectra. 
Due to the general nature of the approach in principle also other mechanisms like DID are contained on a formal level 
in the description, although they cannot be explicitly treated without using a specific theory. 

In our approach we concentrated on a specific mechanism for the coupling of light to the motion of the linear 
molecules, by assuming that the principal axis of the polarizability tensor agree with the principal axis of the inertia 
tensor of the molecule. Under this assumption Eq. (^) and (^are completely general. The importance of the index 
of hclicity m seems to be at variance with the approach of since there only m — components appear. But if 
we express S2\{q,t) by correlation functions of the hydrodynamic variables plus a background spectra with the help 
of the exact generalized constitutive relations as used above we obtain in the q - frame for g ^ 

Sqi{z) = {p2^\R'''"{z)\p2m) 

W{{p2..\R'''"{z)\a..)^f{j\\\R{z)\j\\) (64) 

+2q\{p2^\R'''"{z)\a,^)^f{j^\R{z)\j^) 

kbT 

Here the reduced resolvent R'^^ {z) — Snrn describes the dynamic perpendicular to the hydrodynamic fluctu- 



ations only i.e. the subscript NH indicates, that the projection operator Q, used in Eq. (64) projects on the space 
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of non hydrodynamic variables. Transforming aij to irreducible spherical components, the spectrum (|^) can be 
written 

l''"(z)=T{z)-q'co8\e/2)a'yH{z){3±\R(z)\n) (65) 
where the background spectrum T{z) and the generalized Pockels constant avniz) are given by 

T{z)=g'^{p2,\R'''"{z)\p2o) 

avHiz) ^ g^l^{p2a\R'^"iz)\a2a)^^ (66) 

In leading order in q (i.e. g = 0) we are now allowed to replace the reduced resolvent in Eq. ( |66| ) with the full 
resolvent Similar equations can be derived for the I^^ spectrum. Eq. (|6^) is exactly the form found in [Q. 

The explicit m dependence is not present anymore in ( |65| ) , since we used that at q = correlation functions between 
fluctuations with the same / and m, which do not contain hydrodynamic poles, do not depend on m anymore. But 
it is important to note, that the dynamic coupling to the transversal current fluctuations is not vanishing, only due 
to the existence of an irreducible m = 1 component of the stress tensor. Also in [ p^ this symmetry was implicitly 
used. But after using it, the reduced matrix element {p2i\R'^'^ {z)\a2i) can be replaced by the one for m = at 
g = in Eq. (|6^). The m dependence is replaced by an explicit dependence on the transversal current spectrum, 
which is only present due to the coupling of S22i<},z) to the transversal current fluctuations. We could now further 
aroceed and express the background spectrum and the Pockels constant with the method, explained in the appendices 
by the correlation functions and memory functions, which we used in sections ^ and IV. But since we have 



to evaluate the center of mass current correlation functions and the center of mass density correlation function at 
small but finite values of q (for example to be able to understand the renormalization of the transversal sound 
waves by rotation translation coupling) , we would arrive at exactly the same theory, which we already derived in the 
mentioned sections. The merit of the approach of |Q is to work out the most general form of light scattering spectra 
using only generalized hydrodynamic and generalized constitutive equations for dielectric fluctuations. Thus rigorous 
constraints concerning the appearance of hydrodynamic excitations in different scattering geometries are formulated. 
For explicit considerations of specific scattering mechanisms, as done in this paper, one has in general to include also 
non hydrodynamic variables. 



VI. CONCLUSION 



In this work we developed a microscopic theory of light scattering for linear molecules, concentrating on the direct 
contribution to the spectra. The starting point of our theory is an exact expression for the spectra in terms of 
correlation functions (Eqs. (^ and (||)). It turns out to be important to accurately take the tensorial character of the 
orientational correlation function into account. This is due to the fact that the orientational components of different 
helicity index m transform in general differently under rotations. It is in this context crucial, that the dynamic 
correlation functions S^{q,u;) (contrary to the memory functions) have to be evaluated at small but finite wave 
vectors due to the following reason. The rotational symmetry allows for the dynamic coupling of m = components 
of tensorial densities pim and rotational currents to the longitudinal and the m = 1 components of tensorial densities 
pim and rotational currents to transversal current fluctuations, respectively. Microscopically the coupling is non 
vanishing due to the induction of local stress by the rotation of the molecules. Therefore the hydrodynamic poles 
show up in the corresponding dynamic correlation function S^{q^Lo). Thus we are not allowed to replace them by 
there value at g = 0, where indeed correlation functions for different values of m, but the same I, are equal. Only the 
m = 2 component of S^{q,uj) which does not couple to any hydrodynamic mode and all memory functions, which - 
due to our choice of dynamic variables - do not contain by construction any hydrodynamic pole, can be replaced by 
its value at g = 0. It is the violation of rotational symmetry on the spatial scale of the light scattering experiments 
due to the existence of hydrodynamic modes, which causes the importance of the helicity index m. 

Based on a projection operator formalism, we formulate a microscopic theory for the correlation functions >S';™(g, w) of 
supercooled molecular liquids, which include all possible couplings to hydrodynamic modes. Simplifying the equations 
of the molecular mode coupling theory extended by transverse currents, we demonstrate explicitly, that a qualitative 
description for light scattering spectra near the glass transition can be achieved, which treats correctly the interplay of 
hydrodynamic modes and structural relaxations. We further derive microscopic expressions which give the influence of 
the rotation-translation coupling onto the hydrodynamic poles. It is also shown how other theories of light scattering 
can be expressed by the quantities, which appear in our theory. 
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The equations (Q) and (||) are in principle not restricted to small q-values. They are therefore also valid for the 
interpretation of X-ray spectra, if we give up the restriction to small wave vectors. Since the MMCT, formulated 
in our paper is a theory for all wave vectors, there will be no principal problems to do this. A possible application 



of the restricted theory in [V , would be to compare spectra of different scattering angles and of different scattering 
geometries in order to obtain microscopically relevant quantities like e.g. the rotation-translation coupling. But to 
obtain reliable results, it is of course necessary to give up the Maxwell ansatz and to include (3 relaxation phenomena 
in the memory functions. Further we like to encourage the evaluation of the orientational components for different m 
values from computer simulations. 
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APPENDIX A: COUPLING BETWEEN POLARIZABILITY AND SHEAR STRESS 

To demonstrate how theories of light-scattering, based on various projection operator formalism can in principal 
be expressed by the correlation functions appearing in molecular MCT, we take as an example the set of variables 
from the book of Berne and Pecora and show how their memory functions can be evaluated using the basis set 
of molecular MCT. 

We use the projection operator formalism in Laplace transformed space (see e.g. p^ ). To explain the appearance of 
a Rytov Dip, in the polarizability and one component of the transverse current jx{q^ z) with jx ~ l/-\/2({ji"}oo + 
*{i-i}oo) are chosen as a minimal basis set. z = uj -\- ie \s the complex frequency. We have shown in section 
P that the polarizability for the depolarized light scattering is in the subspace of the 1=2 density. Therefore we 
define avH{Q,z) =: p2 = BmP2m- Applying Mori-Zwanzig projection technique with a projection operator 
Pbp = I1O2) (p^jpa) iP^ \ + I J a:) |j ) {jx\ SLud Qbp =1— Pbp, the dynamics of the polarizability correlation function for 
the depolarized light scattering results from a solution of the following 2x2 matrix equation: 



Z - {p2CR[Cp2) -ip2CR[Cjx) 

~{jxCR[Cp2) z-{jxCR[Cjx 



/ V 



(Al) 



which is still exact. C is the Liouvillian and R'l — Qbp(z~Qbp^Qbp)^^Qbp is the reduced dynamics and (P2P2)'' is 
the static correlation function. The occurrence of a dip in the spectrum relies on the fact that the off diagonal element 
{p2CR'iCjx) does not vanish for small but finite wave-vector q. It follows from momentum conservation, that it is of 
order q. It is therefore possible to define an effective coupling constant R between transversal current fluctuations 
and polarizability fluctuations due to the rotation of the molecule (see ref. |2^ p. 317 to have an explicit connection 
between the phenomenological theory in p^ ] based on an incomplete basis set for the projection technique and our 
microscopic theory). 

R^ lim \im\\{p2CR[C]xf ^0 (A2) 

In the following we show, how this matrix element can be expressed by the memory-functions of MMCT. In this 
theory not only the density p2m, but also the corresponding currents {jq }2m are used as additional variables. A 
coupling of the form Eq. ( [A2| ) can therefore not appear since the reduced resolvent R'2 in the new set of variables is 
projecting perpendicular to the currents i.e also perpendicular to Lp2 as defined above. We have shown in detail in 
sec. 0: 

= ^ B„,Cp2miq, Z)=Y^ B,nq{j^}2rn{q, z) + a/2(2 + l){j^}2m(g, z) 

m 7n 

= J2 B^VW+^Hjohmiq, z) + 0{q) 

m 

=■■ 32 (A3) 
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where T refers to translational currents which occur when applying the LiouviUian on the time dependent positions 
and R refers to rotational currents which appear when applying the LiouviUian to the time dependent orientations. 
For small wave-vectors, we can neglect the contribution of translational currents in the following analysis. We also 
showed in that only the m — 1 component of the rotational current £j2 has non-vanishing matrix elements with 
Cjx- Applying again Mori-Zwanzig projection technique with the enlarged Hilbert space with 



|P2)(P2| + \jx){jx\ + |j2)(j2 



(A4) 



and Ql = X — Pl leaves a 3x3 Matrix equation to be solved. Note that due to Eq. ( |A3| ) contributions that contain 
QlJ~'P2 vanish. 



- ijxCR2Cjx) 
^U^CR'^Cj.,) 



-(P2^J2) 
-(jxCR'2Cj2) 
- {J2CR'2C]2) 





here R'2 is the reduced dynamics due to the new variable set: 

R'^^Qdz-QLCQLy^QL 



(A5) 



(A6) 



Making an additional projection step to obtain a n eff ective 2x2-matrix with \p2) and \jx) as variables (see appendix 
gives a theory for the matrix elements of Eq. (|Al|) : 



{p2CR!^Cp2) = 



\{P2Ci2 



U2CR'2Cj2) 



{jxCR[£jx) = UxCR'2Cjx) 



\{jxCR'2Cj2)\' 

z - {]2CR'2Cn) 



{jxCR[Cp2) 



{jxCR'2Cj2) (j2^P2) 
Z - {j2CR'2Cj2) 



(A7) 



The transformation to the complete set of variables used in MMCT involves inversions of much larger mat rices, but 
the strategy will be the same. For 5 ^ only the matrix elements in Eq. (A7) should be relevant. Eq. ( |A7| ) together 



coefficient Eq. (A2) of the theory described in 



with the mode cou pling expressions in sec. [I], therefore constitute a microscopic theory for the effective coupling 



22| 



APPENDIX B: CONNECTION BETWEEN DIFFERENT PROJECTION SCHEMES 



Lets assume we have a particular basis system A^, Bk with {Bi\Ak) — for all «, /c, which spans the subspace Ti 
An example for that would be the basis set of MMCT described in sec 



Bk the currents {jJJ'}; 



[II, with Ai being the densities pim and 
With such a basis set of the Hilbert space a projection operator can be defined, which 



projects into the subspace 7i^. Within the subspace Ti^ we have a theory to calculate the matrix elements of the 
memory function. In order to be able to compare different projection schemes using a reduced set of variables At, 
which are elements of a subspace TL^ C we need a formalism which expresses all matrix elements in the subspace 
Ti.^ as functions of the matrix elements in the subspace . This can be achieved by applying a formalism described 
in ||]. 

In Ti.^ the operator (2 — £) ^ is given by: 



P^{z- C)P^ - P^CR^'CP'^^ 



(Bl) 



where = 1 — P^ is the usual projector perpendicular to and P^ is the reduced dynamics. 

Due to the orthogonality of the variables Ai, Bk we can decompose P^ in P^ — P^+P^ such that Q^+P^ = 1—P^. 
In the subspace Ti."^ we can write down a similar equation as above: 
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where the reduced dynamics is in the subspace Ti^ is given by 



(B2) 



(B3) 



Let us now express depending on the reduced dynamics R^ . We will therefore have to do an inversion of the 
operator M {Q^iC - z)Q^) 



L' 



M 



pB{C-z)pB pB{L-z)Q^ 
Q''{C-z)pB Q'^{C-z)Q'^ 



The inversion of M gives 



with 



= k'^ ^ R^ - {P^ - R^ M)K{MR^ - P^ 



k = (P^{-z + C)P^ - P^CR^'cP^^ 

Therefore the connection between the two reduced dynamics R^ and R^ are given by: 

R"^' = R^' - (P^ - R.^'cP^)k{P^CR^' - P^) 



(B4) 



(B5) 



(B6) 



(B7) 



This connection was used in Eq. (A7) to derive the connection between the dynamics in the two different basis sets. 



APPENDIX C: MATRIX ELEMENTS OF THE POLARIZABILITY TENSOR 



The polarizability of every molecule Eq. (g) is a tensor of rank 2. It can be written as a scalar plus an irreducible 
tensor of rank 2. In a body fixed coordinate system, with chosen along the principal axis, it has the simple form: 



1 



2ff 



(CI) 



The irreducible spherical components are calculated by: 



a;m = ^ X! C{lll;mim2m)UmriUm23C( 

i,j mi, 1712 



(C2) 



where the i,j € {x,y,z} are Cartesian indices and mi,r7i2 G { — 1,0,1} are spherical "helicity" indices. 
C{lil2l]rnim2m) are the Clebsch Gordan coefhcients. The matrix U is given by 



Ur, 



_!_ n 

1 
-— --^ 



(C3) 



In the body fixed frame only the spherical components with m = do not vanish: = — \/3 a and 
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g. The spherical components in the q - frame are easily obtained by rotation. 



47r 
2/ + 1 



(C4) 



where I?Jj^(rii(t)), Yim{^i{t)) are the Wigner matrices and spherical harmonics, respectively. We used, that the angle 
fli{t) denotes the rotation carrying the body fixed frame into coincidence with the space fixed q - frame. This is the 
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inverse Cti{t) = fl ^{t) to the angle describing the orientation of the molecules with respect t o th e q - frame. The 
Cartesian components in the q - frame are obtained by applying the inverse transformation to (C2) 



"y =E J2 Uj^U^„\C{lll-m^m2m)a 



s 

Im 



(C5) 



The final result for is: 
1 



a^it) = a 



(C6) 



1 



-^Y2o{n{t))+Re{Y22m))) Im(Y22m))) -ReiY^im))) 

Im{Y22{n{t))) ~^Y2omt))-Re{Y22{n{t))) -Im{Y2i{n{t))) 



^ /2vr 

~Re{Y2i{n{t))) -Im{Y2i{n{t))) ^/^Y2o{n{t)) 

where Im, Re denote imaginary and real part, respectively, 

From the matrix element of Eq. (C6) and Eqs. (|^) - (||) we can calculate the contributions which are observable 
in different scattering geometries. Due to the fact that we use the q - frame as the external coordinate system the 
correlation function S*;™™ are diagonal with respect to m,m' and S™,{q,t) — S^,"^{q,t). With 



N 



(o)\ 



J2 {ReY2,n {n, {t))e-'^"' ReY2,n (0))< 



Mrj (0) 



and 



N 



the results (^, (||) are obtained. 
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Laser 




FIG. 1. The possible scattering geometries for a light-scattering experiment. The coordinate system which is denoted with 
{x, y, z) is the common coordinate system used in light scattering theories. The coordinate system which is denoted with 
(a;', y',z') is the q - frame usually used in theoretical descriptions for the dynamics of liquids where the z-axis points along the 
q-vector (arrow) which is probed. 



2.4 




FIG. 2. The depolarization ratio for backscattering geometry D{q,uj) = [{S"l2{<l,'^) + 1/3<S"'22(9, i^)) /S"22i<l,i^) for some 
arbitrary chosen parameters. In units of the rotational frequency uj^ = 1, the other parameters where chosen as cy = 0.6, 
Gs = 1, vr — 5, Ki = 1, Kb ~ 1, KiB^ = Ksr = 1/2. The values of the a-relaxation time r were set identical for all components 

with T = Tl =TR = TSR = TlR = TS e {1, 0.01, 0.001, 10""}. 
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FIG. 3. The depolarization ratio for the same parameters as in Fig. Q on a logarithmic frequency scale. 
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FIG. 5. The susceptibility spectrum x"oo('?i ~ '-^S"oo(q, lo) for the same parameters as in Fig. on a logarithmic frequency 
scale. 



Rj'tov-dipjT = 1 
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FIG. 6. The spectrum 5" 22(13, t^) for the same parameters as in Fig. inset shows the value around zero frequency for 

a = 1 to show the Rytov dip. 
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FIG. 9. The susceptibility spectrum x"22(9>2, 

uj) = uj^^iif,il^RA8 in Fig. I on a logarithmic frequency scale. 
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